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We compute the analogue Hawking radiation for modes which possess a small 
wave vector perpendicular to the horizon. For low frequencies, the resulting mass 
term induces a total reflection. This reflection is accompanied by an extra mode 
mixing which occurs in the supersonic region, and which cancels out the infrared 
divergence of the near horizon spectrum. As a result, the amplitude of the undulation 
(0-frequency wave with macroscopic amplitude) emitted in white hole flows now 
saturates at the linear level, unlike what is found in the massless case. In addition, 
we point out that the mass introduces a new type of undulation which is produced 
in black hole flows, and which is well described in the hydrodynamical regime. 
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Introduction 



Recent studies of the analogue Hawking radiation \T] have shown that a standing (zero- 
frequency) wave is emitted in the supersonic region of white hole flows [21 13] . This wave possesses 
a macroscopic amplitude and a short wavelength fixed by the dispersive properties of the medium. 
Interestingly, it corresponds to some well known solutions in hydrodynamics [4j ^, namely undu- 
lations associated with hydraulic jumps. In addition, undulations have been recently observed in 
water tank experiments [6l [7] aiming to detect the analogue Hawking radiation, but their relation 
with the Hawking effect was not pointed out. This relation was understood in the context of 
atomic Bose-Einstein condensates (BEC), where the emission of an undulation was explained in 
terms of a combination of several effects. Firstly, rico, the spectrum of massless phonons sponta- 
neously produced d la Hawking from the sonic horizon diverges like l/u for — )• 0, as in the Planck 
distribution. Secondly, pu = Pu}=o, the wave number of the undulation is a non-trivial solution of 
the dispersion relation, and thirdly, its group velocity is oriented away from the horizon. 

To see this in more detail first note that the Bogoliubov dispersion relation in a one dimensional 
stationary flow, and for a longitudinal wave vector p, is 



where a; is the conserved frequency, v is the flow velocity, c is the speed of sound, and ^ = ^/2matc 
is the healing length, given in terms of the mass of the atoms mat- The it sign refers to positive 
and negative norm branches, see e.g. [2| for details. The flow profiles giving rise to a black hole 
(BH) and a white hole (WH) sonic horizon are represented in Fig. [TJ The dispersion relation of 
Eq. ([T]) evaluated in the asymptotic supersonic region is plotted in Fig. [2} The zero-frequency roots 
are ±p^, where Pij is 



In this equation, vl and cl are the asymptotic values of the velocity and speed of sound in the 
supersonic region L, and A = I/c^^^l = 2mat/^ characterizes the short distance dispersion. This 
root only exists in a supersonic flow, and its associated group velocity Vgr = l/d^jp is directed 
against the flow. Hence in a BH flow, it is oriented towards the horizon, whereas for a WH one 

^ Even though these are well known, they are not completely understood. In the Introduction of Ref. [5], one 
reads: 'Still, the characteristics and the formation of an undular hydraulic jump are not fully understood'. 



Q = LO — Vp = i: \/^?JP(T+~^JP) , 



(1) 




(2) 
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Figure 1: Shown in this figure are examples of one dimensional black hole flow (solid line) and 
white hole flow (dashed line) with regular asymptotic properties. They are related to each other by 
reversing the sign of the velocity v{x) — )■ —v{x). In both cases, the subsonic R region < c = 1 
is on the right of the horizon, while the supersonic L region is on the left. The near horizon region 
(NHR), where v ^ —1 + kx is a good approximation, has a width in units of k of Dl on the left 
and of Dji on the right. 

it is oriented away from it. This explains why the zero- frequency mode only appears in WH flows 
where it is generated at the sonic horizon. At this point it should be mentioned that these solutions 
are not restricted to superluminal dispersion. A completely similar phenomenon exists in fluids 
characterized by a subluminal dispersion relation, such as that obtained by replacing — t- — ^■^ in 
Eq. Q. This time however, the zero frequency root, and the corresponding undulation, live in the 
subsonic R region of the WH flow. This can be understood because of the (approximate) symmetry 
of the mode equation found in [3j, which replaces a superluminal dispersion by a subluminal one. 

When considering elongated quasi one dimensional systems, but relaxing the assumption that 
the phonon excitations are purely longitudinal, the phonon modes are now characterized by their 
transverse wave number p±, which takes discrete values 2iTn/L±, where n is an integer and L± 
is the characteristic size of the perpendicular dimensions. When 7^ 0, the modified dispersion 
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Figure 2: The solid line represents the positive branch of the comoving frequency Jl, while the 
dashed line represents the negative branch. One clearly sees that the superluminal Bogoliubov 
dispersion is responsible for the two zero- frequency roots ±p^- The sign of the group velocity can 
be seen from the slope of the solid line at the corresponding root. 



relation replacing Eq. ([T]) is 

Q = ^_yp = ±^c^(p2 +pI){i + eip^+pD) ■ (3) 

It is represented in Fig. [sj When p^^^ <^ 1, c^p^ acts as a mass squared. In this regime, there are 
two new zero- frequency roots ±p™. They live in the hydrodynamical regime, characterized by a 
relativistic linear dispersion relation. Indeed, in the limit C'^p'j_ — )• 0, and if v\/(?^ is not too close 
to 1, is independent of ^ and given by 

Pu = ^J^ . (4) 

In addition we note that the group velocity of this new solution has a sign opposite to that of 
Eq. Q. Hence, this new solution will be emitted in BH flows but not in WH ones. 

In brief, we see that the introduction of a perpendicular momentum opens the possibility of 
finding 'massive' undulations in BH flows, which are well-described in the hydrodynamical approx- 
imation of the underlying condensed matter system. To verify if this is the case, one should see how 
the mass affects the spectrum, and in particular if it acts as an infrared regulator that saturates 



5 




Figure 3: As in Fig. [2| the solid line represents the positive branch of the comoving frequency 
while the dashed line represents the negative branch. One sees that p\, which acts as a mass, is 
responsible for new zero-frequency roots ipj^ which occur in the phonon part of the dispersion 
relation. 

the growth of the undulation amplitude found in the massless case [21 [3]. In this paper, these 
issues will be investigated in a simplified context where the phonon modes obey a second order 
differential equation which is a massive Klein-Gordon equation in a curved metric. This analogy 
should work not only for BEC but for other condensed matters systems where the quasi-particle 
dispersion relation is linear at low frequency. We notice that similar issues related to the Cerenkov 
effect have been recently discussed in [8j. 

This paper is organized as follows. In Sec.|lJ we study the solutions of the Klein-Gordon equation 
in a stationary BH metric. We explain how the in /out scattering matrix can be decomposed into 
three blocks that each encode some aspect of mode mixing of massive fields. In Sec. |TIJ we study 
three preparatory cases which are then combined so as to obtain the S'-matrix in a black hole 
flow similar to that represented in Fig. [l} In Sec. |III[ we study the properties of the two-point 
correlation function in the low frequency sector, and its relationship with the undulations. 
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I. SETTINGS, MODE MIXING, STRUCTURE OF THE 5-MATRIX 

In this paper, we study the behavior of a massive scalar field in the 1+1 dimensional metric 

ds^ = c^df - (dx - v{x)dtf. (5) 

For simplicity we impose the condition that the speed of sound is constant, and work in units where 
(? = 1. The stationary metric of Eq. ^ possesses a Killing field Kt = dt, whose norm is It 
is time-like in subsonic flows, space-like in supersonic ones, and it vanishes when 1 — crosses 0. 
At that place, there is a horizon. In what follows, we work with u < and with monotonic flows 
l-v^ that possess a horizon at a; = 0, where v = —1. The definition of the surface gravity of the 
horizon we shall use is 

K=^a,.(l-^;2)|o. (6) 

We adopt this local definition, which no longer refers to the norm of Kt at infinity, because it allows 
us to compare various geometries starting from the near horizon region (NHR). For a black hole, we 
have K > 0, while for a white hole k < 0. Unless specified otherwise, we shall only consider black 
holes. Notice also that we shall work with flow velocities that are either asymptotically bounded 
or unbounded; in the latter case, there will be a singularity. 

The field will be studied at fixed Killing frequency uj = —{Kt)^P^, using a decomposition into 
stationary modes 

^ = j (pMe-'^^'duj. (7) 

At fixed CO, the Klein-Gordon equation gives 

[{uj + id^v){oj + ivd^) + dl- rr?]<P^{x) = 0. (8) 

Similar equations are obtained when studying acoustic perturbations on a fluid flow with a velocity 
profile v{x) |10t I13j. In these cases, a non zero transverse momentum p_L plays the role of the mass 
m. In App. |A| we present the various (effectively) massive wave equations and their differences. 
In what follows, we consider only Eq. ([s]), for profiles v{x) that give rise to analytically soluble 
equations. Yet, we aim to extract generic features. When studying numerically the phonon mode 
equation in a Bose condensate and with a varying sound speed [H] , we recovered the features found 
for solutions to Eq. ([s]). 
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A. Classical trajectories 

To understand the consequences of the mass on black hole radiation, it is useful to first consider 
the corresponding classical problem where p = {dx)^P^ is the momentum of the massive particle 
at fixed oj. In that case, the Hamilton-Jacobi equation associated with Eq. ([s]) is 

^1? = {oj - v{x)pf = p"^ + n? , (9) 

where Q. = u — vp \s the comoving frequency. Eq. ([9|) admits two roots 



p± = -""^V^"'-f (10) 
1 — 



The classical trajectories obey Hamilton's equations dx/dt = Yjd^p and dp/dt = —1/d^x. We 
summarize here their main features with cj > 0, see |15H17j for more details. 



Close to the horizon, at first order in 1 — f ^ ~ 2kx <C 1, one has 



P+ = , (11) 

KX 

= ('2' 

We see that p-^- diverges for x — t- whereas p- hardly varies. The corresponding geodesies 
follow 

x+{t)=x<ie^\ (13) 

x.{t) = 2;° - ^t. (14) 

+ 

The second trajectory is regularly falling across the horizon, while the first undergoes an 
infinite focusing in the past, in a mass independent manner. 

Far away from the horizon, in the left region (L) for 1 — f ^ < 0, both solutions are moving 
to the left (since v < 0) even though < and p_ > 0. 

For 1 — > 0, in the right region (R), as long as (1 — v^)'m? < lo'^, there are two real roots. 
At some point xtp we reach (1 — v'^)m? = where they become complex. This means that 
the trajectory is refiected and falls back across the horizon. Hence, the asymptotic value of 
1 — t;^ > determines the threshold frequency 



ujR = m 



V^-vl, (15) 
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above which the trajectory is not reflected. When u < ojr, there is a single trajectory with 
Pt^ > that starts from the horizon to the right and bounces back across the horizon, see 
FigUl For oj > instead, there are 2 disconnected trajectories, one is moving outwards 
from the horizon, while the other falls in from x = oo. As we shall see, the dimensionality 
of asymptotic modes will be different above and below ojn. 




Figure 4: In this figure trajectories of massive particles for a fixed frequency a; > below the 



threshold of Eq. (15) are shown. The positive momentum trajectory is reflected in the outside 
region at the turning point xtp(ti;) whereas the negative momentum one (or equivalently the positive 
momentum one with negative w [18j) propagates in the inside region where the Killing field is 
space-like. The negative momentum particle has a negative comoving frequency = u — vp, and 
corresponds to a negative norm mode, as we shall see later in the text. 
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B. Mode mixing 

We decompose the field operator in a basis of stationary modes 

^{t, x) = V / ai<liix) + ait(0i^(x))* e-'^^'du; + h.c. , (16) 
J Jo ^ J 

wliere the discrete index j takes into account the dimensionality of mode basis at fixed oj. The 

basis is orthonormal in the sense of the Klein-Gordon scalar product 



dx, 



= ±5{uj - uj')6jj> . (17) 

Following the standard conventions we name the negative norm modes {(p-uj)* so that e*'^*(/)-a; 
is a positive norm mode with negative frequency. 

To obtain the dimensionality of the mode basis, one must identify the solutions of Eq. ([s]) 
that are asymptotically bounded modes (ABM). This requirement univocally picks out a complete 
basis over which the canonical field (j) must be decomposed [25]. Asymptotically, solving the mode 
equation ([s]) is equivalent to solving the Hamilton- Jacobi equation Q. Hence, the dimensionality 
of the ABM can be found by considering the real roots of Eq. ([9]). Moreover, the sign of the norm of 
an asymptotic mode is given by the sign of the corresponding comoving frequency ^}{pi) = oj — vpi, 



as can be seen from Eq. (17). 

In addition, because the situation is non- homogeneous, modes mix and the basis is not unique. 
As usual, we introduce in modes (/)™ and out modes (/>°"* by examining the mode behavior at early 



and late times; see the discussion after Eq. (27) for more precision. The S'-matrix then relates the 
in and out bases. When there is a horizon, these two basis are inequivalent because positive and 
negative norm modes coexist and mix. To further study the mixing, one should consider separately 



the frequencies below and above ur in Eq. (15). 

For < cl! < uJn, there are two ABM. One has a negative norm and propagates behind the 
horizon. The other has a positive norm, comes out from the horizon, and bounces back across the 
horizon, see Fig. |4j The S'-matrix thus has the form 

(18) 

To follow the standard definition of the S-matrix, we use its transpose here. 
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For 00 > upi, there are three ABM. The negative norm one still propagates behind the horizon. 
The second one has a positive norm, comes out from the horizon, and reaches infinity. The third one 
comes from infinity and falls into the hole. As in |25j. we denote the first two with the superscript 
u and the last one with f , because at high momentum, when the mass is negligible, they follow 
retarded (n) and advanced {v) null geodesies. We then define S by 







\ 




/ 


1 out,u 


( 




* 




( 


1 out,u\ 


v 


(puj' 






v 


/ out,v 



(19) 



In this regime, when starting from vacuum, the three out occupation numbers n", and 
obey + = n^u} because of the stationarity of the settings. The first two are given by the 
square of the overlaps 



As we shall see, in both cases, it is useful to decompose the total S-matrix as 

•S* = Sfar • Sext ' 'S'nHRj 



(20) 



(21) 



where each 5- matrix describes one step of the in/out scattering. The first one, S'nHRj describes 
the mode mixing which arises for high momenta p ^ m, near the horizon where the modes are 
effectively massless. The second matrix Sext encodes the elastic scattering which occurs in the 
external region R. Below the threshold, it describes the total reflection, while above it governs the 
grey body factors encoding the partial transmission. The last matrix Sfar describes the mixing 
occurring in the left region between the two modes that are propagating towards x = — oo. 

It should be mentioned that this decomposition is not unique, as only the in/out S-matrix 
is univocally defined. However, in the absence of dispersion, each S-matrix is solution of a well 
defined and independent scattering problem. In addition, this decomposition is very useful as it 
allows us to compute S, and to understand its properties. 



C. Near horizon scattering 



We start with Snhr because its properties are valid for all metrics possessing a horizon and 
because they are determined for momenta much higher than the mass and in the immediate vicinity 
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of the horizon. To simphfy the mode equation (|8|), we introduce the auxihary mode ip^ 



-tu) \ / ,, ax 



Eq. ([8|) is then cast in a canonical form, without the term Hnear in 5; 



-di + 



+ 



We notice that the norm of ip^^ is, up to a sign, given by the Wronksian 



fiuix) = 0. 



(22) 



(23) 



(24) 



Using Eq. (22), one verifies that unit Wronskian 99^ modes, give rise to (j)^; modes that have a 



unit norm with respect to the scalar product of Eq. (17). The relative sign is given by that of the 
comoving frequency Q of Eq. ([9|) . 



In the close vicinity of the horizon, the mass term becomes negligible in Eq. ( 23 ) . More precisely. 



in the near horizon region where 1 — ~ 2kx, keeping only the leading term for kx <C 1, one 



obtains 



1 CO 

- + 



4^2 



(Puj{x) = 0. 



(25) 



Therefore the leading behavior of i^^^ is 



^ui ~ Q{—x)A\2kx\ 

x^O s 



+ e(x) A' \2kx\ 



focusing on the left focusing on tlie riglit 

+ e(-x) B\2kx\-'^^^ + G(x) B'\2kx\' 



(26) 



regularly falling in 



When re-expressing this in terms of the original mode (puj using Eq. ( 22 ) , we see that the B weighted 



terms are regular, and are in fact constant. Therefore, they account for the regularity of the left 
moving mode as it crosses the horizon. Hence, we impose 



B = B'. 



(27) 



On the other hand, the A parts in Eq. (26) oscillate infinitely around x = 0, and account for high 



momentum modes living on either side of the horizon, and which are singular on it. As understood 
by Unruh [20], it is appropriate to combine them in superpositions that are analytic either in the 
upper, or lower, half complex x-plane. The reason is that their analytical properties guarantee 
that these modes correctly characterize the stationary vacuum state which is regular across the 
horizon, and which plays the role of the in vacuum. This characterization applies to the (p^^ modes 
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which are solutions to Eq. ([s]). Hence, the modes (^^j of Eq. (22) are products of a non analytic 
function and an analytic one, which is an Unruh mode: 

\x\-'f^+l X (x + ie)*« , (28) 



r 




e 2k 






V87r2K; 


r 




e 2k 






V87r2K; 



2 X X 



(29) 



We have used Eq. ([24j) to normalize these modes, and their phases have been chosen in order to 
obtain simple expressions. When there are turning points, as is the case for dispersive fields [3] 
and for massive fields, one should pay attention to these phases. 

The normalized modes that propagate on either side of the horizon and vanish on the other side 



are 



i2Kxr2K 



V4 



-Left 



e(-x) 



■KOJ 



|2kx|*2k + 2 



\/47ra; 



The near horizon S-matrix /Snhr is then defined by 

>,NHR 5NHR ' 
?NHR Q,NHR 





(30) 
(31) 

(32) 



Using the analytic properties of the in modes of Eqs. (|28|), (|29|), we immediately obtain 

•Snhr = 



UJTT UJTX 

e2K e 2k 



' e 2k e 2k 



(33) 



Unlike the other matrices in Eq. (21), S'nhr is universal in that it only depends on k. in Eq. (|6j). 
It is independent of the other properties of the profile v(x), and also of the mass ra. In fact, when 
considering a two dimensional massless field, which obeys Eq. (Isl) with m = 0, the left moving v- 



modes decouple, in Eq. (20) vanishes, and the total 5-matrix reduces to the above S'nhr (when 
the asymptotic fiow velocity v is such that out modes are well defined). In that case, on both sides 
of the horizon, the flux of u-quanta is Planckian, and at the standard Hawking temperature k/2'k, 



smce n, 



u 



n"(^ and 



+1 



A 



NHR 



NHR 



-27ra;/K 



(34) 



D. Exterior and interior scatterings 



As mentioned above, 5ext and 5far both depend on other properties of the proflle v{x) than n. 
However their structure can be analyzed in general terms, and the meaning of their coefficients can 
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be identified. Before computing these coefficients in specific flows, it is of value to present their 
general features. 



Below ujfi of Eq. ( 15 ), the positive norm mode is totally reflected while the negative norm mode 



propagates inside the horizon. The exterior scattering matrix S'ext is thus fully characterized by 
the phase accumulated by the positive norm mode in the right region 



r 



(35) 



Since this mode is totally reflected, it is the unique ABM of Eq. (23 ) in R. For small x — )• 0+, using 



Eq. (30), its behavior will be 



~ c 

0+ 



(36) 



|2kx|'2« + 2 + e " X \2kx\ 

which will allow us in the next Sections to extract the phase e*''" . 

In the interior region, there is some extra mode mixing which is described by S'far. Because the 
norms of the two modes are of opposite sign, this mixing introduces new Bogoliubov coefficients: 

(37) 





This scattering is entirely governed by Eq. (23) in the interior region L. When working in an 



appropriate basis, namely when positive and negative frequency modes are complex conjugated. 



the real character of Eq. (23) guarantees that the matrix S'far is an element of SU{\, 1) 



We see that the regularity condition of Eq. (27) for the mode crossing the horizon plus the ABM 



requirement reduces the dimensionality of the four unknown coefficients of Eq. (26) to two. The 



S-matrix in Eq. (|21|) is then 

„Tot oTot ' 



s 



a. 



Tot Q,Tot 




(38) 



This decomposition is depicted in Figjs] Notice that it is the transposed version of the intermediate 



5-matrices of Eqs. (32), (35), (37) that appear in Eq. (38), as explained after Eq. (18) 



For frequencies larger than ujr, the u mode and the v mode, both of positive norm, mix in the 
exterior region R. The matrix Scxt now is 3 x 3 and reads 



iLcft 



r 



10 

Rui Tu) 



i,Leit\ 



,NHR,v 



(39) 
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Figure 5: In this figure a schematic representation of the 5- matrix decomposition in Eq. (38) is 



shown. The asymptotic in and out ampHtudes of the positive frequency mode c/)™ are indicated. 



The 2x2 non trivial sector of this matrix is an element of U{2), and it describes an elastic 
scattering. The interior scattering shares the same properties as occur for frequencies below loh, 
namely the positive and negative norm modes propagating in the region L mix. The outgoing u 
mode is not affected by this scattering, and is thus left unchanged by S'far- Therefore, the structure 
of the total S'-matrix is 

An n \ 



1 

/3^f" 



far ^far 



Tui Ruj 
1 

Ru, T,„ 



\ /^NHR oNHR n\ 



5NHR ;^NHR n 

1 



(40) 



II. EXACTLY SOLVABLE MODELS 



We shall first compute the above 5-matrices in three preparatory cases in order to understand 
various aspects regarding the scattering of massive modes on a stationary horizon. The results 
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of these three cases wih then be put together so as to obtain the S-matrix in a background flow 
relevant for anafogue gravity models, and similar to that presented in Fig. [TJ 

In what follows, the various geometries will be characterized by a single function, given by the 
conformal factor of Eq. ([s]) 

C{x) = l-v^{x). (41) 



The reason to refer only to this function is double. First, as we see from Eq. (23), the mode 
equation for (p^^ only depends on C(x). Second, it will allow us to consider ranges of x where 
C(x) > 1. In such regions, the function v{x) of Eq. ([s]) is complex. However, neither the geometry 



nor the wave equation (23) is ill defined, as can be seen by making the change of time coordinate 



ts = t + J vdx/{l — v'^), which gives 



ds' = il-v')d4-^^. (42) 



This line element depends only on C(x), and can be extended to C > 1. In fact, the v dependent 



prefactor in Eq. (22) accounts for the coordinate change t ^ ts- The status of the auxiliary mode 
iPuj is thus clear: it is the stationary mode when using the 'Schwarzschild' time coordinate ts, which 
is singular on the horizon. While the auxiliary mode Lfuj obeys an equation which is simpler to 



solve, it is singular across the horizon, see Eq. (30). Hence we shall use the original mode (puj to 



impose regularity conditions on the horizon, see Eq. (27). 



Rindler horizon 



It is instructive to first study a Rindler (future) horizon in the above formalism. To do so, we 
use the profile defined by 

C(x) = 2kx. (43) 

It is straightforward to check that this metric has a vanishing scalar curvature R = —d'^C/2, and 
thus describes fiat space. In Eq. (43) k is the 'surface gravity' as defined by Eq. The fact that 
it depends on the arbitrary normalization of the Killing field Kt is free of physical consequence. 



because the S-matrix depends only on the ratio oj/k, see e.g., Eq. (33) 



In this geometry, Eq. (23) reads 



" ^ 2kx 



1 



fLu{x) = 0. 



(44) 
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The interesting aspect of Rindler space is that we know the result in advance. Indeed, since there 



is no pair creation in flat space, the total Bogoliubov transformation of Eq. (21) must be trivial, 
i.e., /33^°* = 0. However, from Eq. (44) we see that close to the horizon the modes behave as in 



Eq. (26), and thus are subjected to the NHR mixing described in Sec IC Therefore, the extra 



scattering described by S'far and 5ext exactly compensates the NHR one, so that the total S-matrix 



is trivial. To show that this is the case, we solve Eq. (44) following the steps of Sec. IC and ID 



Eq. (44) should thus be solved separately for x > and x < 0. On both sides, its solutions can 



be expressed in terms of Bessel functions. We start by studying the exterior R region. For x > 0, 
the only ABM is 



2i 



vr ' \ V 2k I 

where K,y{z) is the Mac-Donald function \21\ and C a constant. At large values of x 

KX ^ ' 



(Puj{x) ~ 2iC . „ 9 9 



e V 2« . 



This exponential decrease is expected since w/j of Eq. ( 15 ) is infinite. Near the horizon, for x 
the ABM behaves as 



0+ 



where 



gfjRindlor X \2kx\ 



— II — -I- — 

2 + 2«;x 2k"^2 



r(l + ia;/K)sinh(^) 



T{iuj/K) ^-2t^ In(f^) 



r(— zw/k) 



(45) 
(46) 

(47) 
(48) 



This is the phase shift that enters in Eq. (35). It will play a crucial role in what follows. 



In the interior region L, the general solution reads 

ipai{x) = A^/^.J_i^/^ 2 



m?x 



2k 



+ B^-xJi^/^ 2\ 



2k 



Near the horizon, for a; — t- , one finds 
^0- ^T{l-iuj/K) 



|2«;xr*2^ + 2 + R 



i2Kxr2«+2 . 



(49) 



(50) 



'2KT{l + iuj/K) 

In order to build the normalized positive frequency left mode (j)]f^^ appearing in Eq. (37), we choose 



A, 



2t\:k 



and = 0, so as to get 



\2kx\ 2k^2 



(51) 



(52) 
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For X — )• — oo, the asymptotic behavior of this mode is |2H I22j. 



Left 



2'KK 



« e 2 X — 



87r2m2 



87r2m2 



(53) 



In the parenthesis, the first term is the asymptotic positive norm out mode, whereas the last factor 



of the second term gives the negative norm one. Therefore, the coefficients of Eq. (37) are 



(54) 

y z, /I r\j 

P^^'* = a^^' X e-f e^^ . (55) 

Making a similar computation for the negative left mode (1;^^^*^^*)*, we obtain a^^" and fi^^'^ and 
verify that Sfar is an element of SU{1, 1). From Eq. ([55]) we see that / a^S'^\'^ = e"^'"^/'^, which 



is exactly the ratio of the near horizon coefficients of Eq. ( 33 ) . This is a necessary condition for 



having = 0. However it is not sufficient, as one also needs the phases to match each other. 



since 



oTot ~NHR«far 



otar I „ tar n 
Pui +<^uj Pi 



far oNHR^j(5i^ 



(56) 



From this equation, one clearly sees the crucial role played by " of Eq. ( 48 ) . An explicit calcu- 



lation shows that the total S'-matrix of Eq. (38) is 

S 



r(— w/k) 







(57) 



e 4 
e"4 

We see that the two in /out coefficients (3^°^ vanish for all values of oo and m. Hence the scattering 
away from the horizon exactly compensates the near horizon mixing and there is no pair creation. 
Of course, this exact cancellation was expected in the present case. However, in more general 
space-times, as we shall see below, a partial cancellation between 5far and S'nhr will be obtained 
for similar reasons. 



B. Totally reflecting model 

Our second example generalizes the former Rindler case in that there is still a total reflection, 
but the profile v{x) now possesses an asymptotically flat interior region. As a result, the asymptotic 
flux of left going particles is well defined, since the emitted quanta are asymptotically described 
by plane waves. The profile which generalizes Eq. ( |43[ ) is 

C{x) = D{-l + e^). (58) 
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The parameter D characterizes the transition from the Near Horizon Region to the asymptotic one. 



In the hmit D — )• oo, C{x) of Eq. (58) becomes that of Eq. (43), which describes Rindler space. 



In the above metric, Eq. (23) is analyticahy solvable in terms of hypergeometric functions [2H 122]. 
The full expression of the general solution is given in App. [Bj To compute the total S-matrix, we 



follow exactly the same procedure as in Sec. II A 



The first important quantity is the phase shift of Eq. (35). To simplify its expression, we 
introduce the dimensionless quantities 



VJ 



CO 

2^' 



2k ' 



The (exact) phase shift then reads 



r(2iro) r(l — ivD + r(l — iw — iQ^ 



„2jroln(D) 



T{-2iw) r(l + m- in+) r(l + iw + 
In the interior L region, the scattering coefficients in 5'far are 

r(l - 2izu)T {-2in+) 



(59) 



(60) 



a, 



far 



2 } Y 



far 



a^^L X e" 



2jroln(D) 



IW 



)r(i 



IW 



im ln{D) 



The total beta coefficient is then given by 



^Tot^~NHR^far + ^far^NHR^..aefl. 



(61) 



(62) 



Its full expression is rather complicated, and not very transparent. It is more interesting to study 
its behavior in different regimes of the parameter space {uj/k, m/ k, D). 

1. Low frequency regime 



An interesting phenomenon happens in the deep infrared regime, a; — )• 0. In this regime, we 



find 



2ft 



2k 



oNHR ^ NHR 



27100 



(63) 
(64) 

(65) 
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These equations show that, while both /3^^^ and diverge as l/w^/^, the total coefficient in 



Eq. (62) does not diverges as l/uj, as one might have expected. Using the analytic character of the 
r functions in Eqs. ( 60|61 ), one finds that the leading term is constant, that is 



(66) 



which is finite for all m > 0. This completely differs from the massless case where /33^°* diverges 
~ P^^^ ~ the discussion after Eq. (33). Moreover, when we take the massless limit of 



Eq. (62) at fixed uj, we obtain the massless result, 

jTot 



f3T ~,/3r", (67) 
for all UJ, and thus in particular we recover the diverging behavior for w — t- 0"*". Before addressing 



the apparent contradiction between Eq. ( 66 ) and Eq. ( 67 ) , it is of value to make a pause and to 



discuss the lesson from Eq. (67). This equation shows that when a massless conformally coupled 
field is scattered on a Killing horizon of a stationary metric which is asymptotically singular in the 
external region (since the curvature -R=— 9^C/2— J-ooforx— )-oo, see Fig. [6]), the particle flux 
is nevertheless well defined in the interior region because it is asymptotically flat, so that the out 



modes of negative Killing frequency are unambiguously defined. In this case, using Eq. (34), one 
gets a Planck spectrum emitted towards asymptotic left infinity. This is rather unusual since the 
Killing frequency is negative, and has in that L region the physical meaning of a momentum since 
the Killing field is space-like. 



The compatibility between Eq. (66) and Eq. (67) is understood when realizing that the saturated 
value Pq"^ of Eq. (66) diverges when m — )■ 0. To see this more precisely, we focus on the regime of 



A 



Tot 



small mass m <C k and small frequencies uj <ti n, for arbitrary ratios uj/m. In this regime, we get 
This expression reveals that there is a change of regime near 



(68) 



ujL = mD2. (69) 

When K » a; » wl, (3"^°^ is growing as in the massless case, whereas for a; ^ it saturates 
at a high but finite value, as can be seen Fig. [Tj As in the case of ultra-violet dispersion [3l [25] . 
we observe that the effective frequency that governs the spectrum depends, as expected, on the 
dispersive frequency, here the mass m, there the ultra-violet scale A, but also depends in a non- 
trivial manner on the parameter D that governs the extension of the NHR. In the present case the 
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Figure 6: Shown in this figure are the Penrose-Carter diagrams of the geometries we shall use. 



On the left, the totally reflecting model of Eq. (58) which is singular in the exterior region, in the 



middle, the CGHS model of Eq. ( 72 ) with its interior singularity, and on the right, the analog model 



of Eq. (82) which is everywhere regular, and can be obtained by pasting the L quadrant of the 



first model with the R quadrant of the second. These diagrams do not represent the full analytic 
extension of each space-time, but only the quadrants that are relevant for our S-matrix, namely, 
the left (L) and Right (R) regions on either side of the (future) horizon %. Precise definitions of 
the various types of infinities along with more details about the last diagram, are given in Ref. |28] . 

power of D is 1/2, whereas for ultra-violet dispersion, the power is (n + l)/n when the dispersion 
relation that replaces Eq. (|9]) is $7^ = +p"+^/A". 

2. Large mass regime 



When the mass is large, i.e., m » k,uj, we find that the coefficients of 5nhr and S'far go to 
their Rindler values in a well controlled manner, e.g., 



gi-jRofl ^ gi<5Rindicr ( I o{K/mD 



This implies that (3'^°^ — )• for m — )• oo as 

P'^^' = 0{K/mDh) = 0{K/uL). 



(70) 



(71) 



This can be understood by considering the Bessel functions of Sec. II A, Their behavior reveals that 



the scattering away from the horizon in region L, which compensates the NHR mixing, occurs on a 
distance from the horizon ~ k/tit?. Hence, for large masses, the entire scattering occurs in a close 
vicinity of the horizon. Therefore, in the large mass limit, the scattering in the present geometry is 
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Tot 1 2 

LO I 




Figure 7: In this figure |/33^°*P is plotted as a function of uo in the regime of low mass and frequencies, 
i.e., uj,m <^ K. For frequencies above the threshold of Eq. (69), behaves as for a massless 



field, and grows as k/u. For frequencies uj < ujl, saturates at ~ k/ujl ~ K/mD^^"^. In the 

opposite regime, where the mass is larger than k/2'k, the suppression arises at a frequency larger 
than the temperature, and 1/33^°* P remains smaller than 1. 

indistinguishable from that occurring in Rindler space. Using WKB techniques [3], which furnish 
reliable approximations in the large mass limit, one can demonstrate that the residual scattering 
outside the NHR is negligible. This means that for m ^ k,u}, irrespectively of the properties of 
the (smooth) profile v{x), the net in/ out Bogoliubov coefficient /^J"* is suppressed by the mass. 



This behavior radically differs from that of the massless case given in Eq. (67), even though both 
cases share the same 5nhr- 



CGHS model 



We now study another exactly soluble example, which is given by the CGHS black hole [9], 
except for the definition of the surface gravity k which is here given by Eq. In Painleve- 
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Gullstrand coordinates, the conformal factor reads 

C{x) = D{1 - e 



(72) 



Even though this geometry is very different from that of Eq. (58) as it is singular in the interior 
region, see Fig. |6j at the level of the mode equation, it gives something very close since the discrete 
interchange C — )• —C and x — t- —x maps one problem into the other. For this reason, the solutions 



of Eq. (23) will also be hypergeometric functions, see App. |B| As in Sec. II B, k is the surface 
gravity and D characterizes the transition from the NHR to the asymptotic region. However, here 



D also controls the value of the threshold frequency lo^i in Eq. (15) since 

ojR = mD^^. (73) 

When ijj < ojR, the positive norm mode is totally reflected, and the accumulated phase shift 
characterizes S'ext- As in the preceding section, to obtain simpler expressions, we introduce 



UJ 



R 



2k 



R 



2k 



which are modified versions of Eq. (59). The exterior phase shift is then 

giicGHs = r(2iro)r {i-iw + n<) r {-iw + n<) ^^^^i^^^^ 

T{-2iw)T [l + iw + Q.<)T [iw + n<) 
From this, conclusions similar to that of Sec. II B can be drawn. For instance, when uj 

recover 

gii^CGHS ^ _\ 



(74) 
(75) 

(76) 

0, we 

(77) 



which is the main ingredient needed to obtain a canceling effect as in Eq. deel), and to have /3,'^°* 



be regular in the limit a; — t- 0. If the mass is large then the behavior is essentially that found for 



Rindler spacetime in Eq. (70), as can be seen by an explicit calculation. 



When UJ > w/j, we are in the configuration where there exist three ABM, as in Eq. (40). 
The greybody factors in the external region R are analytically obtained from the hypergeometric 
functions. The transmission and reflection coefficients are 



W 



3 r f-if]- 



) r (l — iw — iQ 



>) g-iroln(D) 



r(i - 2itu)r (-2fO>) 

r(l + 2iw)T (— zll> — iw^ r (l — iw — ifl>) 
r(l - 2iw)T (-iJ7> + iw) T{l + iw - iCly) ' 

r (2iCi^) r (— iO> — iw) r (i — iw — iCiy) 
r (-2in>) r {in> - iw) r{i-iw + iCiy) ' 



-2iroln(D) 



(78) 
(79) 
(80) 
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Using Eq. ( 40 ) , the asymptotic out-going flux of Eq. ( 20 ) is 



n^ = (0,m|(«r'")^«r'10,m) = |A 



NHR^ 



(81) 



At w = ujR, vanishes, and below lor it is trivially 0. In the next section, the transition shall be 
analyzed in more detail. 

D. Analog model 



We now consider a profile that combines the regular interior region of Sec. II B|with the regular 



C{x) = l-v^{x) 



exterior region of the above CGHS model so as to get a flow similar to that of Fig.[T} The resulting 
geometry is relevant for analog models where the velocity proflle is everywhere bounded. We thus 
consider 

' 2k,x 

Dii-l + e^) for (x < 0), 

^ (82) 
Dr{1 -e for {x > 0). 

This proflle is C^, i.e., it is continuous and its flrst derivative is continuous. This ensures that the 

global geometry obtained is regular, in the sense that the curvature does not contain a distributional 

contribution ^. Since the scattering matrices 5oxt and S'far have been already studied both in the 

exterior and interior regions, all we need to do here is to combine them to get the total S-matrix 



S — ^5fai^ • Sext, • 5nHR • 
SccHTbI SecHTCl SeclTCl 



The two threshold frequencies of Eq. ( 73 ) and Eq. ( 69 ) are now 

1 1 



(83) 



(84) 



ujR = mD^, ujL=mDl. 

Having different values for Dr and will allow us to distinguish their roles. 

We first consider the totally reflecting regime, a; < (jJr. Interestingly, we recover the transition 
seen in Sec. II B and in Fig. [7} Indeed, for uj <^ lor 

gi^cGHS 1 (85) 

Together with the coefficients of S'far, this ensures that has a finite value in the limit a; — t- 0. 
More precisely, in the high k, regime, for m,a; ^ k, we have 

,Tot|2 K,{Dl + Dr) 



(86) 



2'kDr Jijj^ + ujl 



^ A.C. would like to thank Y. Bardoux for fruitful discussions concerning this possibility and for suggesting 
reference [H], which discusses this type of regularity issue. 
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To observe a divergent regime l/S^""*]^ ex k/uj, one needs to assume that col <^ to < lor, where the 
last inequahty is required in order to be below the threshold lor. This implies Dl <^ Dr, hence 

w^'^r - f^^=. (87) 

This expression shows the transition between the diverging regime at the standard temperature, 
which is independent of Dr and Dl, and a saturating regime governed hy ul- For large masses 
^ K,u}, the results are the same as in Sees. II B and II C[ namely the various scattering coefficients 



m 



asymptote to their Rindler values. In numerical simulations [M] all of these results have been 
recovered using rather different settings where the sound speed c varies with x and the velocity v 



is a constant. This demonstrates that the low frequency behavior of Eq. (85) applies to a much 
wider class of situations than the one considered here. 

We now have all the ingredients necessary to study the effects of a massive field on the outgoing 



fluxes when starting from vacuum. On the right side, the outgoing particle flux is that of Sec. II C 
as expected, it vanishes below ujr and above it is given by 

To observe the transition, we work in the high k regime, i.e., u},m <^ k, and obtain 



n^e^—ejuj-LOR) ^ ^ ' ^. (89) 



2tt 



— + u 



The flux is thus continuous when crossing ojr. 

On the left side, the particle fluxes are more complicated since two contributions are present. 



see Eq. (20); is composed of positive frequency particles and nt^ = '^Ji + composed of the 



negative frequency partners. We first notice that both of these are well defined since the profile of 



Eq. (82) is asymptotically flat in L. Using Eq. (40), in full generality, reads 



(90) 

The first term in the first line, which is proportional to Ri^, describes the stimulated production 
in the L region due to the reflected Hawking quanta. The other terms describe the interference 
between the mixing in the NHR region and the scattering in the L region away from the horizon. 
Just as for e^'^'^'^^s below the threshold, the phase of R^j is crucial since there is interference between 
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these two terms. Particularly interesting is behavior of near the threshold frequency loh. In the 
regime of large surface gravity, k ^ u;,m, ioi u > lor > ojl-, one finds 



1 



n,. 



27r 




whereas for lo < uji 



n,. 



K {1 + Dl/Dr) 
2^ 'u^+iol 



(91) 



(92) 



These two equations describe the effects on the spectrum in the L region which are due to a small 



mass. We first notice that is continuous across ujr^ but with a cusp, see Fig. [8j From Eq. (89), 
we see that this is also true for and hence for n!i„, as well. 




Figure 8: Particle flux of positive frequency quanta spontaneously emitted from the horizon towards 
X — )■ — oo in the high k regime. Above wr, one finds a small contribution which is due to the 
reflection (the backscattering) of Hawking quanta emitted towards x = oo. Below the threshold, 
the entire thermal flux is reflected and grows like k/uj for decreasing values of w, until one reaches 



uji where it saturates, as explained in Sec. II B 



We also see that the spectrum depends on the mass m only through the two critical frequencies 



of Eq. (84). It is thus through them that the profile properties, namely and which govern 
the extension of the NHR on the right and on the left, see Fig. [T| affect the spectrum. The lesson 
here is that when dealing with a conformally invariant massless field, the surface gravity is the 
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only background quantity that affects the spectrum. When breaking conformal invariance, by a 
mass, or a non conformal coupling as in 3+1 dimensions, or by adding some ultraviolet dispersion, 
other properties of the background flow affect the spectrum. From the above analysis, and that 
of [3} I25f[27] . the most important ones are the extensions of the near horizon region, on both sides. 

III. MASSIVE MODES AND UNDULATIONS 

In [2], it was noticed that, in a Bose condensate, the density-density correlation function com- 
puted in a stationary white hole flow displays an infrared divergence with a specific short distance 
pattern. As explained in the Introduction, this divergence originates from the combination of three 
properties, firstly the divergence of the Bogoliubov coefficient ~ l/o;^/^ for — )• 0, secondly the 
existence of the zero frequency root p^=o 7^ of Eq. ([2]), and thirdly the group velocity oriented 
away from the sonic horizon. As a result, the zero frequency phonons emitted by the white hole 
give a large contribution to the two-point function which has a short wavelength pattern. More- 
over, this contribution factorizes [3] in the product of twice the same real wave, as in a presence of 
a coherent state of phonon modes . 

In brief, the analysis of the analog Hawking radiation in dispersive media (with either superlu- 
minal or subluminal dispersion [3] ) predicts that white hole flows should emit a zero- frequency wave 
which has a large amplitude and which behaves classically. Interestingly, this type of wave called 
undulation is well known in hydrodynamics [4, 5J. Moreover, undulations have been observed in 
recent experiments aiming to detect the analogue Hawking radiation in water tanks [6l[7j. However 
their intimate relationship to Hawking radiation has not been previously pointed out. Given the 
fact that the low frequency Bogoliubov coefficient /3(^ contributes to the undulation amplitude, we 
think that studying and observing these waves should be conceived as part of the enterprise to 
observe the analogue Hawking radiation. As a last remark, we wish to stress that the linearized 
treatment predicts that, when starting from a vacuum or a thermal state, the undulation ampli- 
tude is described by a Gaussian stochastic ensemble with a vanishing mean amplitude, exactly as 
primordial density fluctuations in the inflationary scenario [30]. Indeed, the mode mixing at the 
horizon amplifies vacuum or thermal fluctuations, and thus generates a random noise. But when 
the initial state is classical, because it contains a large wave packet described by a coherent state, 
the outcome of the amplification is also deterministic. 
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Adding a mass, or a perpendicular momentum, to a two dimensional massless relativistic disper- 
sion relation, also engenders a nontrivial zero frequency root of Eq. ^ in the supersonic L region, 
thereby opening the possibility of having a 'massive undulation'. This possibility has been the main 
motivation of this paper. In the following, we first present the necessary and sufficient conditions 
to find an undulation, starting from the two-point function, as it clearly reveals a subtle aspect, 
namely the cumulative role of low frequencies in determining the amplitude of the zero-frequency 
wave. Then, we study the role of a non zero mass, both in black and white holes fiows. 



A. General properties 



We consider the two-point correlation function evaluated in the in-vacuum 

G{t,t';x,x') = {0,in\(j){t, x)cj}{t' , x')\0,in) . 



Using Eq. (16) we find for equal times that 

G{t,t;x,x') 



Gu}{x, x')duj. 



(93) 



(94) 



In the infrared sector, for frequencies below ujfi in Eq. ( |15[ ), using Eq. (18), the relevant term in 
the integrand is 

G^{x, x') = C(x) (C(^'))* + r Ja:) (r Jx'))* . (95) 
In general, (p"^ and (p™^ are two different functions of x, and G^j is complex and cannot be factorized. 



To see that under some conditions it factorizes, we use Eq. (38) to work with the out modes. Then 



the first term in Eq. ( 95 ) becomes 



C(x) (C(x'))^ 



Tot 1 2 -out 



a,. 



0r (x) (^r (x'))* + 1/3^1' (0°"^(x))* r-^tix) 



+ 2Re {aT°t/3T°t*0-t(x),/.-*(x')} • (96) 
In the limit a; — )• 0, two effects are combined. First (/)°,"* and (/)°"* become the same function of 



(x). This is true in general, but not in the particular case of the two dimensional massless 



field in a black hole metric because in that case (/>°^* and (^°^ vanish on the L and R quadrant 
respectively (see [3]). 



Second, when assuming that > 1 for w ^ 0, since S of Eq. ( |38[ ) is an element of C/(l, 1), 

one has 



Tot|2 



-Tot|2 



5Tot 1 2 



2ie I oTot 1 2 



(97) 
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where e^'^ is a phase. These two facts guarantee that becomes real and factorizes as 

G^{x,x') ~ 8|/3;^°f X cI>u(x)cDu(:e'), (98) 

where the real wave 

$u(x) = Re{e^Vr(^)}, (99) 

gives the profile of the undulation. It should first be noticed that its phase is locked. Indeed, if one 
modifies the arbitrary phase of the out mode (/'q'^*', the modified phase 9 would exactly compensate 
this change so that <I>u would remain unchanged. This can be understood from the fact that <I>u 
oscillates on one side of the horizon and decays on the other. In fact it behaves like an Airy 
function, even though the phase e*^ is different, see [3] for details. One should point out that 
this factorization means that the undulation contributes to observables in the same way that a 
coherent state does, see App. C in Ref. [11]. It thus behaves as a classical wave in that its profile 
and its phase are not random. However, in the present linearized treatment, its amplitude is still 
a random variable, i.e., its mean value is identically zero, and l/?^'"*!^ gives the (w-contribution 
of its) standard deviation. This should be conceived as an important prediction of the linearized 
treatment, and could be validated for a BEC using numerical techniques similar to those of [2], 
and perhaps also in future experiments. 



So far we have worked at fixed u. We now consider the integral over low frequencies in Eq. (94). 
We recall that the divergence of for a; — t- accounts for a growth in time of the r.m.s. 

amplitude of the undulation [21 [3] . When considering an observable evaluated at a finite time t 



after the formation of the horizon, the stationary settings of Eq. (94) with a dense set of frequencies 
should be used with care. Indeed, after such a lapse, one cannot resolve frequencies separated by 
less than 27r/t, as in the Golden Rule. This effectively introduces an infrared cut-off in the integral 
over w. Taking this into account gives the growing rate that depends on the power of the divergence 
of For example, when the Bogoliubov coefficients take their massless values in Eq. (33), in 

the vacuum, the infrared contribution of G grows as 



Gm{t;x,x') ^8 — ^ x <i>^{x) <^uix'), 



4k 

~ — ln(V27r) X ^>u(a;) ^u{x'). (100) 
vr 

Of course, in a medium, this growth would saturate because of the non-linearities, as was observed 
in a Bose condensate [2]. In the experiments of [H [T], only a constant (saturated) amplitude 
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was observed. It would be very interesting to conceive experiments to observe the growth of the 
undulation amplitude. It would be also important to understand if the randomness of the amplitude 
found in the linearized treatment (in vacuum and in thermal states) is replaced by a deterministic 
non-linear behavior when non-linearities are included, or if some residual randomness persists. In 
conclusion, the linearized treatment of perturbations predicts that all (sufficiently regular |26| ) 
white hole flows in dispersive media should emit an undulation with a significant amplitude. 



B. Massive hydrodynamical undulations in black holes 



We now study the undulation in the analog black hole metric of Eq. (82 ). In preceding sections, 
we saw that the low frequency massive modes end up in the inside L region for both signs of u;, 
see Fig. |4j Moreover, in the limit w — )• 0, their momentum (solution of Eq. ([9|, see also Figj3|, is 
finite and given by 

p]^ = p^^, = niD'J (^1 + O (^-^^y (101) 
This means that the zero frequency mode (j)Q^^{x) is a non trivial function of x, opening the 



possibility of finding a behavior similar to that of Eq. (100). When uj^ an explicit calculation 



of Re {e*^(^Q"*(x)}, similar to that made in [3], tells us that the asymptotic profile of $u is 

^^ix) = -^^cos{p^x). (102) 

The next important aspect concerns the calculation of the net contribution of low frequency 
modes to G. In this respect two aspects should be discussed. The first one concerns the fact that 
|^Tot|2 j-^Q longer diverges for w — t- 0. However, the criterion for the factorization of G is only that 



I 1. When ujl k, as shown in Eq. (87), this is the case for frequencies uj <^ k/Itt. The 



second aspect concerns the frequency interval ^ a; < wu such that the momenta are close 
enough to undulation momentum pu so that the out modes t/)"*^* contribute coherently to ^g"*- 
Using Eq. ^ in the asymptotic interior region, we get 

Therefore, at time t after the formation of the horizon, the contribution of the low frequency modes 
to the 2 point function is 

Gin{t;x,x') = 8 / W^°'\^du; x <^^{x) ^^{x'). (104) 

J2lT/t 
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When assuming uu <^ k, which is the case for a small enough mass, using Eq. (87) we get 

''^'^ \P^-'\'du= — \smh-'l , , 



27T/t 27r 



(105) 



Hence, for short times, the amplitude grows as ln(t), as in the massless case. However, when 
t > 2tt/ljl, the amplitude saturates and stays constant afterwards. This is an important prediction 
of this paper. It shows how the transition in oj near ojl with respect to the massless spectrum (see 
Fig. [T]) produces here a change in time of the growth rate. 

To conclude this section, we wish to provide a qualitative evaluation of the importance of this 
infrared contribution to G. To do so, we need to consider some observables, such as the stress 
energy tensor. In particular, its trace accounts for the mass density of the field 



Tr(T) = (r^^)=m^(0^(x)). (106) 
At late times, i.e., t ^ Itt/lol, the contribution of the undulation to the trace is 

TV(fiR) = ^ sinh-i (^) X cos^ ip^ x) . (107) 



Since we work with k ^ lol = mD^^'^, this contribution to the trace is much smaller than k^, the 
typical energy density contained in the Hawking flux for a light fleld. Hence we expect that the 
undulation will not be easily visible in this case. Moreover, when Dl ^ oo, which is the Rindler 
limit, the amplitude goes to as 1/Dl, confirming the stability of Minkowski space. It is also 



interesting to notice that the parameter plays no role (as long as a; ^ un, so that Eq. (87) 
stands), confirming that the undulation is controlled by the interior geometry. 



C. Massive dispersive undulations in white holes 

For white holes, undulations can be found when the dispersion relation is non-relativistic in the 
ultraviolet sector, as discussed in the Introduction. For simplicity, we consider here the superlu- 
minal relation 17^ = + +p'*/A^ which is obtained from Eq. ^ in the limit ^p± — )■ 0. Notice 
however that both superluminal and subluminal dispersion relations give rise to undulations in 
white hole fiows, in virtue of the symmetry which relates them when interchanging at the same 
time the R and L regions, see Sec. III.E in [3j. 
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For superluminal quartic dispersion, the outgoing momentum at zero frequency is found in the 
supersonic region and is given by 



UJ, 



(108) 



(109) 



and the asymptotic behavior of the undulation is [3] 

^A^^^ ^ cos(p^ x + 9u) ^ 

The phase Ojj cannot be obtained from the preceeding equations, because it is mainly determined 
by the dispersive properties of the modes. Using the results of [3], one can establish that, when 
m = Q,9u = (AL»^/^)/(6k) + 7r/4, see also [3l]. 

In the presence of ultraviolet dispersion, the width of frequencies that contribute coherently to 
G is 



(110) 



Since A can be much larger than k, cjy can be either smaller or larger than the Hawking tem- 
perature k/Itt. In what follows we work with oj^ ^ k where the Bogoliubov coefficients are 
well approximated [3] by their relativistic values computed in preceding sections. Therefore, the 
contribution of the low frequency dispersive modes is given by 



GiR(t;x,x') = 8 / \Pl°rdu X $^^(x)$^^(x'), 



27r/t 



(111) 



which is Eq. (104) with and uu replaced by and uj^j. The exact expression of in 



Eq. (62) is quite complicated. To get an undulation, we assume ujl <^ k <^ oj^. In that regime. 



/S^'"* is large for a; < Th, but for a; > Th it becomes exponentially small. Thus, one has 



2lT/t 



lA 



Tot 1 2 



duo 



K./2-K 



2lT/t 



lA 



Tot 1 2 



dijJ. 



In that range of frequencies, /33^°* is well approximated by Eq. (87), therefore 

4k 



Gm.{t;x,x') 



sinh 



K 



sinh ^ ( I 

\0OLtJ 



Hence, at late times and for lo^ S> k, we obtain 

GiB,{t;x,x') = — In 



^'Mx)^Ux'). 



(112) 



(113) 



(114) 



When considering a BEC, the relationship between the scalar field (j) and the density fluctuation 
6p is 5p oc dx4> [29]. Hence the mean value of the equal-time density-density two point function is 



{d,^{x)d,,^{x')) = ^ln 



nuJL 



sin (^pijx + djj) sin (pf/x' -|- 9, 



U 



(115) 
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This generalizes what was found in [2] in that, in the supersonic region, one stih finds a short 
distance checker board pattern in the x, x' plane, and the undulation amplitude still grows initially 
as ln(t). However, when there is a mass term, it grows only for a finite amount of time ~ 2Tr /ujl, 
after which it saturates. The mass therefore provides a saturation mechanism that can occur before 
nonlinearities take place. Moreover, because oc A S> k, the r.m.s. amplitude of the undulation 
is large. 

So far we have considered only the case where the initial state is a vacuum. When dealing with 
a thermal state, as discussed in [2j the initial growth rate is no longer logarithmic but linear in 
time. However, the mass term acts again as an infrared regulator because the initial distribution 
of phonons is expressed in terms of J7 > m, and not in terms of the constant frequency cj 
Hence no divergence is found when integrating over to when computing the two point function. In 
addition, the random character of the undulation amplitude is fully preserved when taking into 
account some initial thermal noise. What is modified is the undulation r.m.s. amplitude. When 
the initial temperature Tin is much larger than Tjj = k/Itt, the above two-point function is, roughly 
speaking, multiplied by T\^JTh- 



IV. CONCLUSIONS 

In this paper we have studied the consequences of a mass term in the mode equation on (the 
analogue of) Hawking radiation. We showed that the scattering of massive modes on a stationary 
black hole horizon is rather complicated. It contains several mode mixing terms which occur on 
various length scales and which interfere with each other in a non-trivial way. In what follows we 
summarize our main results. 

Notwithstanding the fact that the mass does not affect the near horizon mode mixing, as can 



be seen in Eq. (33), the mass does regularize the infrared behavior of the net in /out Bogoliubov 



transformation in Eq. (83). The reason for this is the extra mode mixing, described by 5'far of 
Eq. (37), which occurs in the supersonic inside region, and which interferes with the near horizon 
scattering so as to cancel out the divergence in of the coefficient. Indeed, the squared 

norm of the total Bogoliubov coefficient saturates as |/33'°*P ~ KjlTiijJi^ for a; — t- 0, where ujl is the 



threshold frequency of Eq. (84). As a consequence, the undulation r.m.s. amplitudes now saturate 



after a lapse of time ~ 2tt/ujl and then stay constant. This has to be contrasted with the massless 
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case where the saturation of the ampHtude can only occur because of non-Unearities, or dissipation, 
in the system. 

The presence of a mass induces a new type of undulation in the supersonic region that exists 
in black hole flows. Unlike the undulations occurring in white hole flows which are due to some 
ultraviolet dispersion, this new type occurs in the hydrodynamical regime if the mass term is small 
enough. It will thus appear both in superluminal and subluminal media. However, as shown by 



Eq. (107), the typical energy density carried by an undulation is small, and thus this new type 
should be difficult to detect. 

Although the S'-matrix coefficients governing black hole and white hole flows are the same, in 
virtue of the transformation v — >• —v that maps one case onto the other, the frequency ranges 
that contribute to the massive and the dispersive undulations are very different as can be seen 



by comparing Eq. (103) with Eq. (110). As a result, the white hole undulations possess larger 
amplitudes. In addition, since the wave length of the undulation is smaller in the white hole case. 



it gives rise to even larger amplitudes for a BEC, as can be seen from Eq. (115). These results 



might also be relevant for surface waves where 'transversal instabilities' have been observed |33| . 
Finally, we emphasize that the properties of the spectrum and the undulations depend on the 



mass m essentially through the effective frequencies ujl and coji of Eq. (84). These frequencies 
are both proportional to m but also depend in a non trivial way on Dl and which determine 
the spatial extension of the near horizon region, on the inside and on the outside respectively. 
Therefore, as in the case of ultraviolet dispersion, these two quantities should be conceived of as 
the most relevant geometrical properties, after the surface gravity k. 
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Appendix A: General acoustic d'Alembert equation 



In the hydrodynamical approximation, i.e., when neglecting short distance dispersion, hnear 
density perturbations in a moving fluid obey a AD d'Alembert equation, in a curved space-time [H 



--d„ 



0. 



The effective metric is 

„2 _ P r„2 ,.2 



\(?dt^ — {dx — Vxdt)^ — {dy — Vydt)^ — {dx — f^dt)^] , 



(Al) 



(A2) 



where p is the density of the fluid, c the sound speed and v the velocity of the background flow. 

When assuming that the flow profile is one dimensional and homogeneous in the perpendicular 
dimensions, the mode equation becomes effectively two dimensional because the transverse wave 
number is constant. In what follows, its norm is noted p±. Assuming in addition that the fiow is 
stationary, at fixed frequency uj, the mode equation reads 

2 



-((J + idxv)^{uj + ivdx) + -dxpdx - pi 
p p 



ix) = 0. 



(A3) 



In the body of the paper, we assumed p{x) = c(x) = 1. This violates the continuity equation of 
fluid mechanics ^ but it allows us to study exactly soluble models. 



It is also interesting to see how Eq. (23) is altered by varying c(x) and p{x). Modifying the field 



redefinition of Eq. (22) 



<pLo{x) 



/X v(^x' )dx' 



p{c'^—v^) 
^2 



-Vuj{x), 



we obtain 



'2 c2 



y2 |^g2 _ ^2^2 



ifujix) = 0, 



where the effective potential is 



Vcix) 



7i 



7i 



(A4) 



(A5) 



(A6) 



We see that this equation has exactly the same structure as the simplified equation (23) that was 
used for some of the calculations in this paper. 



^ Note that this condition does not concern surface waves because the effective metric is different from that 
given above [71 [32] . 
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Appendix B: Hypergeometric solutions 



We provide the exact solutions of Eq. ( 23 ) for the profiles considered in Sec. II B and Sec. II C In 



both cases, the solutions are hypergeometric functions. For definiteness, we recall their definition 

^-^ i (aji (oji (c + nj n! 



where a, h and z are complex numbers, and c is a non negative integer. Their main properties, and 
the asymptotic behaviors of these functions can be found in [" 21} [22] . 



We start with the profile of Eq. (58), and use the notation of Eq. (59). For x > 0, the only 



ABM is 

(^^(x) = C(e^ - 1)^ (1 - e-^y^e~^F {l - iw + in+, 1 - iw - 2; e'^rr ) , 
where C is an arbitrary constant. For x < 0, the general solution is 

(Puj{x) =A{1 - e^)^e~^ - (iw - in+,iw + in+; 1 + 2iw; 1 - e"^ 



+ B{1 - e^)^(e"^ - (-iw - iQ+, -iw + in+; 1 - 2iw; 1 - e"^ 



with A and B arbitrary constants. 



For the profile in Eq. (72), and for x > 0, the general solution is 



2kx , 1 , 2kx 



(Puj{x) =A{1 - e^~)2{e^ - ly^F \^iw - iil,iw + iVl;l + 2iw; 1 - 

+ B{1 - e"^)^(e^ - (-iw - iU, -iw + iO; 1 - 2iw; 1 - e"s" 

Here A and B are arbitrary constants and = for u > ul and = iO< for lo < lol. The 



definitions of these dimensionless quantities are given in Sec. II C 
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